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Data Loading Method, Transmitter, and Base Station 
Field 

The invention relates to a data loading method in a communication 
system, a transmitter and a base station. 

5 Background 

Systems based on using multiple antennas are one of the most 
promising techniques for achieving high data rates. By combining MIMO (mul- 
tiple input multiple output) techniques with OFDM (orthogonal frequency divi- 
sion multiplexing) modulation, the frequency selective MIMO channel is turned 
10 into a set of frequency flat MIMO fading channels which can be individually 
processed. 

If a channel is known at a transmitter side, the channel matrix can 
be decomposed for each sub-carrier using SVD (singular value decomposi- 
tion). As a result a set of orthogonal sub-channels is obtained in the space 
15 domain. These elementary sub-channels are also called eigenmodes. 

One prior art method to find the optimal bit and power allocation for 
a set of multiple parallel sub-channels (for example, for eigenmodes of a 
MIMO-OFDM system) is Hughes-Hartog (HH) algorithm. Hughes-Hartog algo- 
rithm is depicted in John A. C: Bingham: Multicarrier Modulation for Data 
20 Transmission: An Idea Whose Time Has Come, IEEE Communications Maga- 
zine, pp. 5-14, May 1990, which is taken herein as a reference. However, the 
required computational effort of Hughes-Hartog algorithm increases with the 
average bit rate and therefore it is not a practical solution for high bit rate sys- 
tems 

25 Some sub-optimal fast loading algorithms have been proposed in 

Peter S. Chow, John M. Cioffi and John A. C. Bingham: A Practical Discrete 
Multitone Transceiver Loading Algorithm over Spectral Shaped Channels, 
IEEE Trans, on Communications, vol. 43, no. 2/3/4, pp. 773-775, Febru- 
ary/March/April 1995, Robert F. H. Fischer and Johannes B. Huber: A New 

30 Loading Algorithm for Discrete Multitone Transmission, IEEE Proceeding of 
Global Telecommunications Conferences, vol. 1, pp. 724-728, Nov 1996 and 
Yu Wei and John M. Cioffi: On Constant Power Water-Filling", IEEE Proceed- 
ing of International Conference on Communications, vol. 6, pp. 1665-1669, 
June 2001 which are taken herein as a reference. 



These prior art methods are developed for slow varying channels 
like asymmetric digital subscriber line. For such channels, the signalling over- 
head required, for instance, for informing a receiver about the modulation 
scheme of each eigenmode can be neglected due to possibility to update the 
transmission parameters at relatively long time intervals. However, when the 
channel is changing fast, as is the case in wireless systems, the transmission 
parameters can be kept constant only for a time period shorter than the chan- 
nel coherence time and therefore the amount of signalling overhead is signifi- 
cant. 

Brief description of the invention 

According to an aspect of the invention, there is provided a data 
loading method in a communication system where sub-carriers include eigen- 
modes, comprising: estimating a channel matrix; calculating a singular value 
decomposition of the estimated channel matrix for obtaining eigenvalue esti- 
mates; defining biases between eigenvalues and eigenvalue estimates and 
performing a channel estimation reliability test based on the defined biases; 
carrying out bias compensation for eigenvalue estimates based on the defined 
biases; calculating equivalent power gain; arranging eigenmodes into a prede- 
termined number of clusters, each cluster comprising eigenmodes of different 
quality levels; pre-allocating transmission power to the eigenmodes according 
to their capacity by using the calculated equivalent power gain; determining 
collective transmission power to be allocated to each cluster based on the pre- 
allocation, selecting the optimum modulation and coding scheme and allocat- 
ing collective transmission power to the eigenmodes. 

According to another aspect of the invention, there is provided a 
transmitter of a communication system where sub-carriers are divided into ei- 
genmodes, comprising: means for estimating a channel matrix; means for cal- 
culating a singular value decomposition of the estimated channel matrix for 
obtaining eigenvalue estimates; means for defining biases between eigenval- 
ues and eigenvalue estimates and performing a channel estimation reliability 
test based on the defined biases; means for carrying out bias compensation for 
eigenvalue estimates based on the defined biases; means for calculating 
equivalent power gain; means for arranging eigenmodes into a predetermined 
number of clusters, each cluster comprising eigenmodes of different quality 
levels; means for pre-allocating transmission power to the eigenmodes accord- 
ing to their capacity by using the calculated equivalent power gain; means for 



determining collective transmission power to be allocated to each cluster 
based on the pre-allocation; means for selecting the optimum modulation and 
coding scheme and allocating collective transmission power to the eigen- 
modes. 

According to another aspect of the invention, there is provided a 
transmitter of a communication system where sub-carriers are divided into ei- 
genmodes, configured to: estimate a channel matrix; calculate a singular value 
decomposition of the estimated channel matrix for obtaining eigenvalue esti- 
mates; define biases between eigenvalues and eigenvalue estimates and per- 
forming a channel estimation reliability test based on the defined biases; carry 
out bias compensation for eigenvalue estimates based on the defined biases; 
calculate equivalent power gain; arrange eigenmodes into a predetermined 
number of clusters, each cluster comprising eigenmodes of different quality 
levels; pre-allocate transmission power to the eigenmodes according to their 
capacity by using the calculated equivalent power gain; determine collective 
transmission power to be allocated to each cluster based on the pre-allocation; 
select the optimum modulation and coding scheme and allocating collective 
transmission power to the eigenmodes. 

According to another aspect of the invention, there is provided a 
base station of a communication system where sub-carriers are divided into 
eigenmodes, comprising: means for estimating a channel matrix; means for 
calculating a singular value decomposition of the estimated channel matrix for 
obtaining eigenvalue estimates; means for defining biases between eigenval- 
ues and eigenvalue estimates and performing a channel estimation reliability 
test based on the defined biases; means for carrying out bias compensation for 
eigenvalue estimates based on the defined biases; means for calculating 
equivalent power gain; means for arranging eigenmodes into a predetermined 
number of clusters, each cluster comprising eigenmodes of different quality 
levels; means for pre-allocating transmission power to the eigenmodes accord- 
ing to their capacity by using the calculated equivalent power gain; means for 
determining collective transmission power to be allocated to each cluster 
based on the pre-allocation; means for selecting the optimum modulation and 
coding scheme and allocating collective transmission power to the eigen- 
modes. 

According to another aspect of the invention, there is provided a 
base station of a communication system, where sub-carriers are divided into 



4 



eigenmodes, configured to: estimate a channel matrix; calculate a singular 
value decomposition of the estimated channel matrix for obtaining eigenvalue 
estimates; define biases between eigenvalues and eigenvalue estimates and 
performing a channel estimation reliability test based on the defined biases; 
5 carry out bias compensation for eigenvalue estimates based on the defined 
biases; calculate equivalent power gain; arrange eigenmodes into a predeter- 
mined number of clusters, each cluster comprising eigenmodes of different 
quality levels; pre-allocate transmission power to the eigenmodes according to 
their capacity by using the calculated equivalent power gain; determine collec- 
10 tive transmission power to be allocated to each cluster based on the pre- 
allocation; select the optimum modulation and coding scheme and allocating 
collective transmission power to the eigenmodes. 

Embodiments of the invention are described in the dependent 

claims. 

1 5 The method and system of the invention provide several advan- 

tages. In a preferred embodiment of the invention, a sub-optimal low- 
complexity bit and power loading method that requires low signalling overhead 
is provided. Therefore the algorithm is especially suitable for high bit rate sys- 
tems. The algorithm is also robust against channel estimation errors at the 

20 transmitter side. 

List of drawings 

In the following, the invention will be described in greater detail with 
reference to the preferred embodiments and the accompanying drawings, in 
which 

25 Figure 1 shows an example of a communication system, 

Figure 2 is a flow chart, and 

Figure 3 illustrates an example of a transmitter structure. 

Description of embodiments 

With reference to Figure 1, we examine an example of a data 
30 transmission system in which the preferred embodiments of the invention can 
be applied. The present invention can be applied in various wireless communi- 
cation systems based on MIMO-OFDM. One example of such a communica- 
tion system is IEEE 802.11 a wireless LAN communication system. The basic 
idea of OFDM is to split a high-rate data stream into parallel streams that are 
35 transmitted simultaneously over different orthogonal sub-carriers. An OFDM 
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signal consists of a sum of sub-carriers that are modulated by using phase 
shift keying (PSK) or quadrature amplitude modulation (QAM). MIMO systems 
include multiple transmission and reception antennas. It is clear to a person 
skilled in the art that the method according to the invention can be applied to 
5 systems utilizing different modulation methods or air interface standards. 

In MIMO-OFDM systems, the frequency selective MIMO channel is 
turned into a set of frequency flat MIMO channels which can be individually 
processed. This reduces computational complexity. Elementary orthogonal 
sub-channels of each sub-carrier, obtained by using SVD (singular value de- 
10 composition) of a channel matrix at each sub-carrier, can also be called ei- 
genmodes. 

The embodiments are not, however, restricted to the system given 
as an example but a person skilled in the art may apply the solution in other 
systems provided with the necessary properties. 

15 Figure 1 is a simplified illustration of a digital data transmission sys- 

tem to which the solution according to the invention is applicable. This is a part 
of a cellular radio system, which comprises a base station or an equivalent 
network element 100, which has bi-directional radio links 102 and 104 to sub- 
scriber terminals 106 and 108. The subscriber terminals may be fixed, vehicle- 

20 mounted or portable. The base station includes transmitters, for instance. From 
the transceivers of the base station there is a connection to an antenna unit, 
which establishes the bi-directional radio links to the subscriber terminal. The 
base station is further connected to a base station controller or an equivalent 
network element 110, which transmits the connections of the terminals to the 

25 other parts of the network. The base station controller controls in a centralized 
manner several base stations connected to it. 

The cellular radio system can also communicate with other networks 
such as a public switched telephone network or the Internet. 

If a transmitter has some knowledge about CSI (channel state in- 

30 formation) the spectral efficiency can be increased by applying adaptive 
transmission techniques. In a TDD (time division duplex) based system, the 
channel can be estimated at the base station on the basis of a signal received 
in one or more selected previous uplink timeslots. The estimate may be used 
for link adaptation during the next down link timeslot. The link adaptation con- 

35 sists of modifying transmission parameters according to channel variation in 
order to maximise the throughput using at most the maximum transmission 



6 



power and fulfilling requirements set for the reliability of the transmission. Typi- 
cally, the reliability is evaluated in terms of frame error rate (FER) or bit error 
rate (BER). 

The transmitter typically informs the receiver about the chosen 
transmission parameters, for instance selected eigenmodes, constellations, 
powers allocated to the eigenmodes and channel code parameters, using a 
signalling channel. The signalling overhead reduces effective capacity of the 
downlink channel, thus the target is to minimize it for making the usage of 
available capacity more efficient. 

Next, by the aid of Figure 2, let us examine in further detail an em- 
bodiment of the data loading method (or algorithm). The example system used 
here includes two receiver antennas and at least two transmitting antennas. 
The scheme exploits radio channel reciprocity of TDD (time division duplex) 
transmission and operates in an open loop adaptation mode. Typically, for 
each transmitted frame, in order to maintain the selected frame error rate 
(FER) under the total transmitted power constraint, the modulation scheme 
and the allocated power to each eigenmode are adapted according to channel 
conditions. 

If a channel is known at a transmitter side, the channel matrix can 
be decomposed for each sub-carrier using SVD (singular value decomposi- 
tion). As a result a set of orthogonal sub-channels is obtained in the space 
domain. These elementary sub-channels are also called in this application ei- 
genmodes. 

For such a scenario, the optimal signalling in terms of channel ca- 
pacity is usually eigenmode transmission, where a complex symbol is allocated 
to each eigenmode and the transmission power is allocated to the eigenmodes 
according to, for instance, a prior art water filling algorithm or modified water 
filling algorithm. The modified water filling algorithm is explained later in further 
detail. The process of adapting power and the amount of information to be al- 
located to each eigenmode is called here bit and power loading. The amount 
of information to be allocated to the eigenmodes varies according to the modu- 
lation level such as complex signal constellation and/or the channel coding 
parameters such as code rate. Ml MO-OFDM eigenmode transmission is de- 
scribed in Kai-Kit Wong, Roger S. K. Cheng, Khaled Ben Letaief and Ross D. 
Murch: Adaptive Spatial-Subcarrier Trellis Coded MQAM and Power Optimiza- 
tion for OFDM transmissions, Proceeding of IEEE Vehicular Technology Con- 



ference, vol. 3, pp. 2049-2053, 2000, H. Sampath, P. Stoica and Arogyaswami 
J. Paulraj: Generalized linear precoder and decoder design for MIMO channels 
using the weighted MMSE criterion, IEEE Trans on Communications, vol. 49, 
no. 12, pp. 2198 -2206, December 2001. 

The embodiment of the method provides reduced signalling over- 
head in relation to the number of sub-carriers by grouping eigenmodes into 
preferably, two clusters and using the same modulation and coding scheme 
(MCS) for all selected eigenmodes belonging to the same cluster. The clusters 
preferably include the strongest and the weakest eigenmodes of each sub- 
carrier. 

The embodiment starts in block 200. 

In block 202, a channel matrix is estimated. 

An estimate of the channel matrix (H c ) at each sub-carrier c is pro- 
vided: 

H C =H C+ N C (1) 
wherein 

c means a sub-carrier and 

H c means the estimated channel matrix and 

H c means the true value of channel matrix and 

N c means the means channel estimation error with the entries 

modelled as independent Gaussian random variables with zero mean and 
variance a 2 N . 

Channel estimation is typically made according to selected prior art 
method. The channel estimation is usually based on pilot signals of the re- 
versed frame. Variance of the error of the channel estimation (a 2 N ) is typically 

also estimated using an appropriate prior art method. Channel estimation is 
well known in prior art and therefore it is not explained here in further detail. 

In block 204, singular value decomposition (SVD) of the estimated 
channel matrix is calculated. As a result the eigenvalues of the estimated 
channel matrix are obtained. Singular value decomposition is calculated typi- 
cally independently for each sub-carrier. 

The SVD of the estimated channel matrix is given by: 
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H c = U A V , 



(2) 



wherein 

c means sub-carrier and 
5 U c is an RXR (R = the number of receiver antennas equal with 2 in 

this example) unitary matrix contains in its columns left singular vectors of the 
channel matrix and 

V c is an TXT (T = the number of transmitter antennas > 2 in this ex- 
ample) unitary matrix contains in its columns right singular vectors of the chan- 
10 nel matrix and 

V C H means the transpose and complex conjugate (hermitian) of the 
matrix V c and 

A c is an RXT (R = the number of receiver antennas equal with 2 in 
this example, T = the number of transmitter antennas) dimensional matrix hav- 
15 ing the elements ^/^7 on the main diagonal and all other elements being 
zeros; the elements on the main diagonal are descending ordered 

The matrix V c is also used for linear pre-filtering of the signal 
20 transmitted at a sub-carrier c in order to obtain eigemodes. Due to the channel 
estimation errors orthogonality between the resulted eigenmodes in each sub- 
carrier is usually lost. Therefore a spatial equalizer is required at the receiver 
side for each sub-carrier. The spatial equalizer can be implemented either in 
zero forcing or linear minimum mean square structure. The equalization is well 
25 known in prior art and therefore it is not explained here in further detail. 

If the channel estimation is fairly accurate (channel estimation errors 
smaller than 10%), the instantaneous signal to noise ratio at the equalizer out- 
put can be approximated by: 

30 SNR j>c «^^ c |v» c v, >c | 2 , (3) 



wherein 

c means sub-carrier and 
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i means the spatial sub-channel (eigenmode) index for a given sub- 
carrier and 

N 0 means the power spectral density of noise at the receiver side 

and 

P i c means the power allocated at the eigenmode (*,c) and 

A i c means the square of the / 'th singular value of the true channel 

matrix H c . 

v i c is the i"th right singular vector of the estimated channel matrix 
H c and is given by the /'th column of matrix V c and 

v" means the transpose and complex conjugate (hermitian) of the 
vector v i c and 

\ i c is the / 'th right singular vector of the true channel matrix H c . 

v" v f c | in equation (3) represents the loss in SNR due to 

the channel estimation errors. It's instantaneous values are difficult to calculate 
but it was found by simulations that: 

E ".jw ! i)- «> 

wherein 

c means sub-carrier and 

i means the spatial sub-channel (eigenmode) index for a given sub- 
carrier and 

e n c { } means the statistical expectation operator, where the expec- 
tation is taken with respect of probability density function of N c elements (N c is 
the channel estimation error) and 

X ic means the square of the /'th singular value of the true channel 
matrix H c and 

X ic means the square of the / 'th singular value of the estimated 
channel matrix H„ . 



Based on relations (3) and (4), the average SNR (signal-to-noise ra- 
tio) loss can be compensated by increasing the instantaneous allocated power 
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X 

to each eigenmode with the ratio . To implement this we need a reliable 

estimate of the \ c . Eigenvalues are squares of the singular values of the 
channel matrix obtained by SVD. 

In practise, the channel matrix itself cannot determine. Instead of 
the real channel matrix, the estimate is used. Correspondingly, also for eigen- 
values estimates are used. 

The main idea of the method is to partially compensate the bias be- 
tween i, c and Aj c . Notice that A i c are the eigenvalues of the hermitian matrix 
H C H" and \ c are the eigenvalues of the hermitian matrix H C H" . 



Let's denote the entries in the 2X2 hermitian matrix H H H as fol- 



lows: 



H H H = 



(5) 



wherein 

a o b c>K> c c represent matrix elements and 
* means a complex conjugate. 

By solving the characteristic equation of the matrix H C H" , the ei- 
genvalues are given by: 

Kc = \ («c + c c + J(a c -c c ) 2 + 4\b c \ 2 \ 

1 / / 5 n ( 6 ) 

wherein 

a c ,b c ,c c are defined in (5) and 
I I means absolute value. 



The estimated channel matrix H e multiplied with its' transpose 
complex-conjugate (Hermitian) matrix is expressed as follows: 
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H H H = 

c c 



a c b c 
b\ c. 



(7) 



wherein 

a c ,b c ,b c ,c c represent matrix elements and 
* means a complex conjugate. 

The eigenvalues of H C H" are given by: 



%c 2 



1 



+ 4 



\c=^\o c +c c -^j(a c -c c ) +4\b c 



(8) 



wherein 

a c ,b c ,c c are defined in (7) and 
I I means absolute value. 

In block 206, biases between eigenvalues l (C and their estimates 
i, c are defined for channel estimation reliability test and the channel estima- 
tion reliability test is performed. 

The biases between eigenvalues X i c and their estimates A i c , the 
terms of expressions (6) and (8), can be calculated by: 



B,=E K {a c -a c } = Tcr 2 N 
B 2 =E Nc {c c -c c }=Ta 2 N 

5 3=E Nc {|4f-|6 c | 2 } = ^|H c t + r^ 

* 4 =E Nc j(« c -c c ) 2 -{a c -c c ) 2 } = 4a* |H e £ + 2To* 



(9) 



wherein 

B X ,B 2 ,B 3 ,B 4 are the bias terms and 
a c ,b c ,c c are defined in (5) and 
a c ,b c ,c c are defined in (7) and 
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E N { } means the statistical expectation operator, where the expec- 
tation is taken with respect of probability density function of N c (N c is the 
channel estimation error) and 

T is the number of transmission antennas and 

<y 2 N is variance of the error of the channel estimation and 

|H e ||p means the square of the Frobenius norm of channel matrix 
H c and it is given by ||H e |£ = a c +c c . 

From (9) is clear that for computing bias terms B 3 and B 4 we need 
first to estimate the square of the Frobenius norm of the channel matrix H c . An 
instantaneous consistent estimate is given by: 



K c \\ ? =a c +c c -B x -B 2 =a c +c c -2T(Tl (10) 

Based on this estimate a test for the channel estimation reliability 
evaluation is carried out for each sub-carrier as follows: 



*(Wf>0 



(11) 

the sub-carrier c will be used in the loading algorithm 
else 

the sub-carrier c will be discarded from the loading algorithm. 

In block 208, bias compensation is carried out for eigenvalue esti- 
mates X, c . 

For the eigenvalue estimation the bias terms given by (9) are used 
to partially compensate the bias in the terms of expression (8) as follows: 



^-2 



1 



= \[{o c - b x y + (c c -B 2 y- p c -c c y + 4\b c \ 2 ^ 



(12) 



wherein 

B } ,B 2 are the bias terms given by (9) and 
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(x) + means the positive value of x defined as: (x) + = max(x,0) and 
B 3 and B 4 represents the estimates of B 3 and B A ; they are ob- 
tained by replacing ||H C |£ in (9) with its instantaneous estimate ||H C ||J given by 
(10). 

For smaller eigenvalues (l, c in equation (12)), the terms, which are 
under the square root are left uncompensated. This produces a slight underes- 
timation of the eigenvalue causing increasing in system robustness. The ei- 
genmodes corresponding to the eigenvalues which become negative after bias 
compensation, are discarded from the loading algorithm. 

In block 210, the equivalent power gain is calculated. 

By using the estimated eigenvalues (12) in (4) and combining with 
(3) we can express the equivalent power gain at the eigenmode (i,c) as fol- 
lows: 

X 2 

S',c=i^ ' (13) 

Clearly, in case of perfect channel estimation at the transmitter side 
the equivalent power gain become the square of the singular value of the 
channel matrix (g ic -——^ j as is known from the prior art. 

An eigenmode is an elementary sub-channel created within one 
sub-carrier by using SVD of the channel matrix at that sub-carrier. The eigen- 
modes are obtained by pre-filtering the transmitted signal with the right hand 
singular vector matrix. 

The power required at the eigenmode (i,c) in order to obtain a 
given target signal to noise ratio, SNR, , is given by: 

/> C (SNR,) = SNR,A (14) 

In the following, we will describe the bit and power loading part of 
the data loading method, using as an example a system with two receiving an- 
tennas and at least two transmitter antennas. The number of eigenmodes ob- 
tained in each sub-carrier is given by the minimum between the number of 
transmit and receive antennas (two in the considered example). 
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Briefly, the main idea of the method is to exploit the intrinsic spatial 
diversity achieved by using MIMO-OFDM transmission. In the following we 
present one embodiment based on turbo channel coding but it is clear to a 
person skilled in the art that turbo codes may be replaced by other prior art 
channel coding methods. 

In block 212, eigenmodes are arranged into a predetermined num- 
ber of clusters, each cluster comprising eigenmodes of different quality levels. 

In one embodiment, 2C eigenmodes (C meaning the total number of 
sub-carriers) are grouped into two clusters, consisting of the strongest and the 
weakest eigenmodes from each sub-carrier, respectively. Except some eigen- 
modes corresponding to the most faded sub-carriers, a small difference be- 
tween the eigenmodes gains belonging to the same cluster is experienced due 
to spatial diversity. Consequently, by skipping some of the most faded eigen- 
modes, the same MCS (modulation and coding scheme) can be used for all 
selected eigenmodes belonging to the same cluster and the required signalling 
overhead is C times reduced. 

In one embodiment, each cluster is independently encoded by using 
single input single output (SISO) turbo code. The resulted bits are usually also 
interleaved and modulated. For each cluster the encoding is performed jointly 
in time and frequency domain, a codeword covering the selected eigenmodes 
from one cluster during whole transmitted frame. In this way the codeword is 
enlarged to achieve interleaving gain while the adaptability between clusters is 
still preserved. 

The modulation and coding parameters (for example the complex 
constellation used in modulation, the coding rate) are modified according to 
channel conditions in order to maintain the target FER (frame error rate). Spe- 
cifically, there are available a set of K distinct modulation and coding schemes 
(MCS), and we denote them by b k , k = \,...,K, the spectral efficiency being 
provided by k'th MCS. 

For each available MCS, the required SNR to achieve the target 
FER in additive white Gaussian noise (AWGN) channel can be obtained ana- 
lytically or by preliminary simulations and the results can be stored into a look 
up table. 

Let's denote by y k , k = \,...,K the SNR required by the rth MCS to 
achieve the target FER in AWGN channel. We assume that b k are given in the 
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ascending order 6, <b 2 <••• <b K and for an efficient MCS set is required that 

Y\ <r 2 < --<Yk- 

For each transmitted frame, first, the total power is divided between 
the clusters according with their capacities, and then, the MCS and the se- 
5 lected eigenmodes in each cluster are independently optimized to maximize 
the throughput subject to equal SNR constraint for all selected eigenmodes. 

Thus, in block 214, transmission power is pre-allocated to the ei- 
genmodes according to their capacity by using the calculated equivalent power 
gain. 

10 In practise, the amount of information allocated to an eigenmode is 

limited not only by the noise power but also by the maximum rate achieved by 
the most spectral efficient MCS (modulation and coding scheme). Therefore, in 
power allocation, instead of prior art water filling a modified water filling method 
can be used for taking into account the saturation of signal-to-noise ratio of the 

15 eigenmode and the average signal-to-noise ratio gap. 

In the modified water filling, the power to be pre-allocated for an ei- 
genmode is expressed: 



20 



Z^ MWF) =min 





r n ] 


X' 













(15) 



wherein 

MWF means modified water filling 

i is a cluster index, 

c is a carrier index, 
25 N 0 means spectral density of noise at the receiver, 

G is the average signal-to-noise ratio gap between the Shannon ca- 
pacity and the spectral efficiency provided by the MCS set, 

(x) + means the positive value of jc defined as: =max(jc,0) and 

y K means the SNR required by the most spectral efficient MCS to 
30 achieve the target FER in AWGN (Additive White Gaussian Noise) channel 

g ic means the equivalent power gain at the eigenmode (i,c) given 

by (13) 

// means the "water level" and it is found from the maximum trans- 
mit power constraint: 
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1=1 c=l 



(16) 



wherein 

P T means the maximum transmit power 
C means the total number of sub-carriers. 

When the modified water filling method described above is used, 



less power is allocated to some eigenmodes than if the prior art water filling 
method was used. These eigenmodes are called saturated eigenmodes be- 
cause although the amount of allocated power can be increased, the amount 
of information transmitted using these eigenmodes has already reached the 
maximum e.g. is saturated. 



In block 216, based on the power pre-allocation of block 214, collec- 



tive power to be allocated to each cluster is determined. The collective power 
allocated to a cluster can be expressed as follows: 



wherein 

C means the number of sub-carriers, 

c is the current sub-carrier, 

2 means a summing operation and 
^(mwf) js obtained from equation (15). 

In block 218, transmission power is allocated to the eigenmodes 



and the optimum modulation and coding scheme (MCS) is selected for each 
cluster. 



The clusters are independently optimised by searching for the MCS 



that provides maximum throughput when the collective power P, is allocated to 
the eigenmodes according to the target FER (frame error rate), under equal 
signal-to-noise constraint. This is, in practise, typically equivalent with the 
maximisation over the MCS set of the product between the maximum number 




(17) 
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of eigenmodes selectable using a certain MCS and the spectral efficiency pro- 
vided by that MCS. 

The maximisation is preferably applied in the same manner to each 
cluster. Therefore the cluster index (/) is omitted for brevity. By denoting with P 
the collective power allocated to the cluster, the index of the optimum MCS, k 0 , 
which maximises the throughput under the total transmitted power and the 
maximum FER constraints can be written as: 



( 



k 0 = arg max b k 



max \S\ 

ceS J 



= arg max b k s k 

k 



(18) 



wherein 

S is a subset of the set of eigenmode indexes with cardinality |S|, 
s k represents the maximum number of eigenmodes that can be se- 
lected using the k'th modulation scheme, 

P c {n) is given by (14) (y k 

C means the number of carriers, 
c is the current carrier, 

b k means the spectral efficiency (number of bits per transmitted 
complex symbol) provided by k'th MCS. 

The values s k can be calculated using the following algorithm: 

Stepl: arranging in a descending order the equivalent eigenmode 
power gains g c . The arranging process means mathematically that the indexes 

vector 

i = [JiJ2*—Jc] • (19) 
such that g Jt >g h >---> g k is found. 

Step2: computing the vector w =[w 15 w 2 ,...,w c ], with each entry given 



by 



w -=^oZ— - (20) 



„=> g Jm 
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wherein 

N 0 means the power spectral density of noise at the receiver side 

and 

c is the index of the entry and , 
n means the summation index, 

j„ means the n'th entry in the indexes vector j given by (19) 
g Jm means the equivalent power gain of the j n 'th eigenmode in the 
current cluster 

£ means the summing operation. 

Step 3: Let s 0 =C and^ o =0. For k = \ \oK, let s h =s k _ x and 
^ w sJ k >P. decrease the number of selected eigenmodes, s k =s k -\, until the 
power constraint, w s j k <P , is satisfied. Every time when the power constraint 
is satisfied the corresponding value s k is saved. 



Using the values s k provided by the previously presented algorithm, 
the optimum MCS , k 0 , is found by using (18) and the powers allocated to the 
cluster's eigemodes are computed based on (14) as follows: 



p jSrO = r^ if n<s ka 

0 if n > j. 



(21) 



wherein 

j'„ means the n'th entry in the indexes vector j given by (19) 

P Jm means the power allocated to the j n 'th eigenmode in the current 

cluster 

k 0 means the index of the selected MCS in the current cluster 
y K means the SNR required by the selected MCS (k 0 ) at the cur- 
rent cluster to achieve the target FER in AWGN channel 

\ means the number of the selected eigenmodes in the current 

cluster 

If correlation between antennas is strong and/or signal-to-noise- 
ration is low, it is possible that the number of selected eigenmodes at the 
weakest cluster is small and the resulted code word is very short. In such a 
case, the performance of codes may decrease and the target FER (frame error 
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rate) is not maintained. In this case, the weakest cluster is not used and the 
power pre-allocated to the weakest cluster is reused by the strongest cluster. 
The arrow 224 depicts this procedure. 

The method ends in block 220. The arrow 222 depicts one possibil- 
ity of repeating the method. 

Next, an example of a transmitter structure according to an em- 
bodiment of the data loading method described above is depicted in further 
detail by the aid of Figure 3. It is obvious for a skilled person that the structure 
of a transmitter may vary from what is depicted in Figure 3. 

The present invention can be applied in various wireless communi- 
cation systems based on Ml MO-OFDM. One example of such a communica- 
tion system is IEEE 802.11 a wireless LAN communication system. In Figure 3, 
there is depicted only a part of an OFDM transmitter focusing on the transmit- 
ter blocks required by the loading method. It is obvious to a person skilled in 
the art that the transmitter may also include elements other than those illus- 
trated in Figure 3. 

The transmitter described here is thought to be constructed of two 
parts: of pre-processor 318 and of OFDM-modulator 320. The data loading 
method described above with the aid of Figure 2 is carried out mainly in bit and 
power loading block 304 which is a part of the pre-processor 318. 

The embodiment of the data loading method described above is 
based on that the transmitter has prior "knowledge of the radio channel. The 
radio channel can be estimated at a receiver being coupled to the transmitter 
on the basis of one or more signals received in selected previous timeslots. 
The link adaptation consists of modifying transmission parameters according to 
channel variation in order to maximise the throughput using at most the maxi- 
mum transmission power and fulfilling requirements set for the reliability of the 
transmission. Typically, the reliability is evaluated in terms of frame error rate 
(FER) or bit error rate (BER). 

The receiver being coupled to the transmitter is not depicted in Fig- 
ure 3 for the sake of clarity. The receiver is in this embodiment a prior art 
MIMO-OFDM receiver. 

The transmitter typically informs the receiver about the chosen 
transmission parameters, for instance selected eigenmodes, constellations, 
powers allocated to the eigenmodes and channel code parameters, using a 
signalling channel. 
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Channel estimation block 300 carries out channel estimation for sin- 
gular value decomposition SVD which in turn is carried out in block 302. Both 
the channel estimation and SVD are described above with the aid of Figure 2. 

Encoding and interleaving blocks 306A-306B carry out channel cod- 
ing, interleaving and modulation. In one embodiment, channel coding is carried 
out by encoding each cluster independently by using a single input single out- 
put (SISO) turbo code. For each cluster the encoding is performed jointly in 
time and frequency domain, a codeword covering the selected eigenmodes 
from one cluster during whole transmitted frame. In this way the codeword is 
enlarged to achieve interleaving gain while the adaptability between clusters is 
still preserved. The resulted bits are usually also interleaved and modulated 
using some prior art interleaving and modulation methods. 

In the linear pre-combining blocks 310A-310B, the complex symbols 
generated in the blocks 306A-306B are scaled according to the powers allo- 
cated to each eigenmodes in block 304 and then the transmitted signal at each 
sub-carrier is filtered with the right hand singular matrix V c (see Figure 2, block 

204). The number of pre-combining blocks is equal to the number of sub- 
carriers. 

In OFDM-modulator of the example of Figure 3, an IFFT is carried 
out in IFFT blocks 312A-312B. It is obvious to a skilled person that .the number 
of blocks may vary according to the implementation. Typically, it is equal to the 
number of transmit antennas. Then, in blocks 314A-314B, the signal is con- 
verted from parallel to serial form and a cyclic prefix is added. 

The OFDM-signal is then conveyed via a DAC (digital-to-analogue 
converter) to the antenna 316A-316B. The number of antennas vary according 
to the implementation. In MIMO systems, there are several antennas. 

OFDM-systems are well known in the art. 

Even though the invention is described above with reference to an 
example according to the accompanying drawings, it is clear that the invention 
is not restricted thereto but it can be modified in several ways within the scope 
of the appended claim. 



